Band Structure and Optical Absorption in Multilayer Armchair Graphene 
Nanoribbons: A Pariser-Parr-Pople Model Study 
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I. INTRODUCTION 



Kondayya Gundra^'0 and Alok Shukla^ 
I ^Department of Physics, Indian Institute of Technology, Bombay, Mumbai 400076 INDIA 

, Using the tight binding and Pariser-Parr-Pople (PPP) model Hamiltonians, we study the elec- 

{^JQ' tronic structure and optical response of multilayer armchair graphene nanoribbons (AGNRs), both 

^ ' with and without a gate bias. In particular, the influence of the number of layers (n), and the 

, strength of the electric field applied perpendicular to layers, for different types of edge alignments, 

■ is explored on their electro-optical properties. As a function of increasing n, the energy gap initially 

CO ' decreases, eventually saturating for large n. The intensity of the linear optical absorption in these 

systems also increases with increasing n, and depends crucially on the polarization direction of the 
incident light, and the type of the edge alignment. This provides an efficient way of determining the 
nature of the edge alignment, and 7i, in the experiments. In the presence of a gate bias, the intensity 
O ' of optical absorption behaves in a nontrivial way. The absorption becomes more intense for the large 

fields in narrow ribbons exhibiting a red shift of the band gap with the increasing field strength, 
while in broad ribbons exhibiting a blue shift, the absorption becomes weaker. However, for smaller 
electric fields, the absorption intensity exhibits more complicated behavior with respect to the field 
strength. Thus, the effect of the gate bias on optical absorption intensity in multilayer AGNRs is in 
sharp contrast to the bilayer graphene, which exhibits only enhancement of the absorption intensity 
"j^ ' with the increasing electric field. 

PACS numbers: 78.20.Bh, 78.67.Wj, 73.22.Pr , 78.40.Ri 

o 

Since the isolation of graphenei, rapid advances have been made in the experimental^ and theoretical"^ inves- 
C\| ] tigations of this truly two dimensional (2D) material, and related nanostructures, for use in the next generation 
^ • opto-electronic devicesi^ii While graphene is a zero gap semiconductor, energy gap opens up in quasi one-dimensional 
t ■ (ID) graphene nanoribbons (GNRs) due to reduced dimensions. The mechanism of gap opening in GNRs depends 
on the nature of the edge termination. First principles calculations^!^ show that, in case of mono-layer zigzag GNRs 
(ZGNRs) energy gaps open up due to edge magnetism. On the other hand, the gap opening in mono-layer armchair 
GNRs (AGNRs), is believed to be both due to quantum confinement, and the reduction of the bond lengths at the 
^ ■ edges ,^ 

While numerous theoretical studies of the electronic structure and related properties of the mono-layer GNRs 
existfiS relatively fewer calculations on bilayer and multilayer-GNRs have been performed i^i"— Recently, gated bilayer 
graphene has attracted a great deal of attention in the experimental;^"— as well as theoretical communities;^"— The 
energy gap in bilayer graphene, in the presence of a transverse electric field, has been found to be tunable over a wide 
range of values (up to 250 meV)^'2^ Furthermore, the electrical noise levels in bilayer graphene channels is much less 
^ , compared to the noise levels in single layer graphene,— thus making bilayer graphene, a promising candidate for the 
fabrication of high-quality electronic devices. Based upon this, one can argue that compared to mono-layer GNRs, 
bilayer GNRs will possess superior properties from the point of view of device physics. Moreover, it is interesting to 
study multilayer graphene and GNRs in their own right, as they can help us in understanding the evolution of the 
electronic structure from graphene to bulk graphitcj^ 

Recently, we have developed an approach to study the electronic structure of graphene nanostructures, based 
upon the 7r-electron Pariser-Parr-Pople (PPP) model Hamiltonian;^ and used it to study the band structure, edge 
magnetism, and optical absorption of mono -layer GNRs of various types, at the Hartree-Fock leveli^^ The advantage 
of the PPP model based methodology is that it incorporates the long-range Coulomb interactions among the tt- 
electrons in a natural way. In this work, we apply our approach^^ to perform a detailed investigation of the electronic 
structure of multilayer AGNRs, with and without an external electric field, and for various edge alignments. For 
the gated ribbons with the intrinsic band gaps above a critical value (say, Cc), the gap decreases with the increasing 
field strength, while for those with gaps lower than ec, it increases with the external field, a result in good agreement 
with the recent ab initio results of Sahu et al^ In addition to the studies of the band structure, we also present 
calculations of the optical absorption spectra of various multilayer AGNRs for different edge alignments, for various 
polarization directions of the incident light, and with, and without, a gate bias. To the best of our knowledge, no prior 
calculations of optical absorption in multilayer-GNRs exist, which take electron-electron interactions into account. 
Our calculations reveal that for the gated bilayer ribbons, for large bias fields, the optical absorption intensity increases 
for the ribbons with the intrinsic band gaps higher than ec, while for the ribbons with the gaps smaller than ec, the 



absorption intensity decreases. However, for smaller bias fields, the absorption intensity exhibits a more complicated 
behavior with respect to the field strength, a behavior in sharp contrast with the gated bilayer graphene which exhibits 
only increase in the absorption intensity with the increasing bias fieldj^'^>^ Furthermore, we find that in addition 
to the polarization direction, the absorption intensity is found to be crucially dependent on the number of layers in 
the multilayer AGNRs, and their edge alignment, thereby allowing for optical determination of their structure. 

The remainder of this paper is organized as follows. In section HTJ we briefly describe our PPP model based 
theoretical methodology. Next, in Sec. IIIII we present the results of our calculations on the band structure of various 
multilayer AGNRs, and discuss the variation of the band gap with the increasing number of layers, their widths, and 
as a function of the gate bias. In section ITVl we present the optical absorption spectra of these systems, and discuss 
their variation with n, external electric field, and the type of edge alignment. Finally, in Sec|V]we summarize our 
results, and present the conclusions. 



II. THEORETICAL METHODOLOGY 



In our earlier works we have used the PPP model to extensively to study the electronic structure and optical 
properties of finite 7r-electron systems such as conjugated molecules and oligomers at various levels of theory^^ while 
in our recent work we extended it to study GNRs in the infinite length limit, with ID periodic boundary conditions 
imposed)^ The PPP model Hamiltoniauf^E with one 7r-electron per carbon atom, is given by 

H = - Ujictcja + 4<tC»<t) + 

i i<j 

where cl^ creates an electron of spin a on the orbital of carbon atom i, Uicr = c\^Cia is the number of electrons with 
the spin a, and rii = riic is the total number of electrons on atom i. The parameters U and Vij are the on-site 
and long-range Coulomb interactions, respectively, while tij is the one-electron hopping matrix element. On setting 
U = Vij = 0, the Hamiltonian reduces to tight binding (TB) model, while on setting just Vij — 0, it reduces to the 
Hubbard model. The parametrization of Coulomb interactions is Ohno like^ 

where, k^j- depicts the dielectric constant of the system which can simulate the effects of screening, and Ri j is the 
distance in A between the i-th and the j-th carbon atoms. In our earlier work on GNRs;^ we used the ab-initio 
GW band structure of mono layer AGNR-12 (AGNR-TVa, denotes an AGNR with Na dimer lines across the width) 
reported by Son et al.~ to obtain a set of "modified screened Coulomb parameters," with U = 6.0 eV and Kij = 2.0 
{i 7^ j) and Ki,i ~ 1. Note that our modified screened parameters are slightly different from the screened parameters 
reported initially by Chandross and Mazumdar^ with U = 8.0 eV and Kij = 2.0 (i 7^ j) and K^.i = 1, aimed 
at describing the optical properties of phenyl-based polymers within the PPP model. With our modified screened 
parameters, we obtained excellent agreement between the ab initio GW band gaps;^ and our PPP Hartree-Fock band 
gaps for mono layer AGNRs of various widths;^ therefore, we have used these parameters in the present study of 
multilayer AGNRs as well. 

In these calculations, we consider Bernal packed multilayer AGNRs with ABAB... repeat pattern, and two possible 
edge alignments, called a and /3 alignmentSfH shown in Fig. [1] The intra layer nearest-neighbor (NN) distance was 
taken to be 1.42 A, and, for the inter layer separation the value 3.35 A, identical to that in graphite, was used. As far 
as the hoppings are concerned, the intra layer NN hopping was chosen to be i = 2.7 eV, along with the next-nearest- 
neighbor (NNN) hopping t = 0.27 eV. For the inter-layer hopping also, we considered not only NN hopping t± — 0.4 
eV, but also the NNN hopping integral tj_ = 0.3 eVi^ We further explore the parameter dependence of the calculated 
optical absorption spectra in section HVBI Because multilayer AGNRs are closed-shell insulating systems, like their 
mono layer counterparts, we used the restricted HF (RHF) method to obtain the band structure in our calculations, 
as described in our previous worki^ 
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Figure 1: The structures of bilayer AGNR (a) a-alignment and (b) /3-alignment 
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Figure 2: (Color online) Variation of the energy gap with the increasing width for bilayer AGNRs, obtained by the tight-binding 
method for (a) a-alignment (b) /3-alignment. 

III. ENERGY GAPS AND BAND STRUCTURE 
A. Energy gaps 

The TB calculations on mono-layer AGNRs^ predicted that AGNR-iV^ with Na = 3p and iV^ = 3p + 1 (p 
is a positive integer) exhibit energy gaps (E^*), which are inversely proportional to Na, whereas AGNRs with 
Na = 3p + 2 are gapless. Hence, based upon TB theory, they can be classified into three groups with energy gaps 
varying as E^p > > EgP~^'^{= 0). However, ab initio density- functional theory (DFT) calculations^ on these 

ribbons predicted a different relationship Eg^'^^ > E^p > E^p^'^{-=/^ 0), with the important result that even for 
Na = 3p + 2, AGNRs exhibit nonzero energy gaps, due to the fact that the bond lengths involving the edge atoms 
are shorter than those in the interior. When the decrease in the bond length is incorporated in the TB approach by 
increasing the corresponding hopping, one also obtains finite gaps for Na = 3p + 2 ribbons, although the relation 
E^P > E^P+^ > E'^P^'^ still holds. Based upon first principles DFT calculations Sahu et al.r^^ discovered that for 
bilayer AGNRs, and other multilayer AGNRs E^/+'^ > E^p > E^p+'^{^ 0) is valid, in full agreement with the DFT 
results of Son et al^ obtained for the mono-layer AGNRs. In what follows we first explore the same relationship using 
the TB model, and in Figs. [2] and [H] we present the variation of the energy gaps of different families of bilayer and 
eight layer AGNRs, respectively, with respect to their widths, both for the a and the /3 alignments. The values of of 
the hopping integrals used in these calculations were as specified in section [Hi except that the hoppings at the edges 
were increased by 12%^ leading to finite band gaps for ribbons with Na = 3p + 2 for the narrow multilayer ribbons. 
As is obvious from the figure, our results are consistent with the aforesaid relation E^p > EgP~^^ > E^p^"^ valid for 
mono-layer AGNRs at the TB level. Furthermore, for a given width of the bilayer and eight layer AGNRs. band gaps 
with the /? alignment are always more than that for the a alignment. 

Next, for each of the three families of AGNRs, we examine the variation of the gap as a function of the increasing 
number of layers (n) at the TB level, and results of our calculations are presented in Fig. 2^ for the a alignment, and 
in Fig. 5)3 for the /3 alignment. Although, these calculations were performed for p = 3, a similar behavior was also 
observed for other values of p. We find that: (a) again the band gaps for the ribbons in /3-alignment are larger than 
those for the a-alignment, and (b) band gaps for both types of alignments saturate fairly rapidly with respect to n, 
suggesting that a small number of layers are needed to achieve the bulk values. 

In particular, for the a-aligned ribbons with Na = 11 (3p -I- 2 class) , the band gap decreases rapidly with n and 
becomes negligibly small for large values of n (Figdlja), blue/dashed line) whereas, for the /3 alignment the band gap 
is much less sensitive to n and attains a constant value, for smaller values of n (FigHl^b), blue/dashed line). This 
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Figure 3: (Color online) Variation of the energy gap with the increasing width for the AGNRs with eight layers, obtained by 
the tight-binding method for (a) a-alignment (b) /3-alignment. 
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Figure 4: (Color online) Variation of energy gap with number of layers obtained by the tight-binding method for multilayer 
AGNRs in (a) Q-alignment (b) /^-alignment. 



tendency holds even when the e-e interactions are considered, as is evident in the Fig. [5] which presents the variation 
of band gap with number of layers for multilayer AGNRs with A^^ = 8 in a-alignment (black line) and in /^-alignment 
(red line), obtained by PPP-RHF approach. From the results it is obvious that the inclusion of the e-e interactions 
leads to the widening up of the band gaps, compared to the TB approach. 



B. Band structure 



For the sake of brevity, in what follows, we denote a given multilayer AGNR as n-AGNR-A^'^-a (n-AGNR-A'^i-/?), 
implying an AGNR with n layers, each of width Na, arranged in a (/3) alignment. In our previous work we have 
emphasized the role of long-range e-e interactions in widening the band gaos of the mono-layer AGNRs)^ For example, 




Figure 5: (Color online) Variation of energy gap with number of layers obtained using the PPP model, within the RHF 
approximation, for the multilayer AGNRs with Na = 8, in a and /3 alignments. 
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Figure 6: (Color online) Band structure of (a) 2-AGNR-ll-a , (b) 2-AGNR-11-/3, obtained using the PPP-RHF approach. 
The inset contains the magnified band structure near the Fermi energy {Ef = 0) of: (a) 2-AGNR-ll (black lines) and (b) 
mono-layer AGNR-ff (red lines). 

we demonstrated that for mono-layer AGNR-11, the Hubbard model based calculations predict a negligible gap, while 
the PPP-RHF method predicts a gapped system in agreement with the GW calculations We expect the long-range 
e-e interactions, as incorporated in the PPP model, to play a similar role for multilayer AGNRs as well. According to 
the ab initio DFT calculations of Sahu et alr^J^, all n-AGNR-iV4-a, of the family JV4 = 3p + 2, exhibit much smaller 
gaps compared to the other families, suggesting a metal-like behavior. In Figs. [6^ and [7] we present the PPP-RHF 
band structures of a few members of this family, 2-AGNR-II-q;, 3-AGNR-8-Q;, 4-AGNR-8-Q;, and lO-AGNR-8-a, and 
from the figures it is obvious that all these ribbons have substantial band gaps. Therefore, we believe that the ab 
initio results of Sahu et alM^ reporting these systems as almost metallic, are due to the well-known tendency of the 
DFT to underestimate the gaps, and that a better estimate of the gap can only be made by some electron-correlated 
approach such as the GW approximation^ 

We investigate the effect of edge alignment on the band structure of bilayer AGNRs by presenting the band structure 
of 2-AGNR-ll-Q! (Figini (a)), and 2-AGNR-11-/3 (FiglHlJb)). The band structure near the Fermi energy is magnified 
in the insets of those figures, which also contain the band structure of mono-layer AGNR-11. In the presence of the 
second layer, each mono-layer band is split into two bands: one with the lower energy and the other with the higher 
energy, thus reducing the gap. From insets of Figs. and[6]3, it is clear that the valence and the conduction bands 
formed due to the aforesaid band splitting for 2-AGNR-II-q; near Ep are separated by larger energy compared to 
the case of 2-AGNR- 11-/3. Therefore, it suggests that the perturbation introduced by the second layer is larger in 
the case of the a-aligned AGNRs, as compared to the /?-aligned AGNRs, for which we offer the following geometrical 
explanation. Intuitively speaking, the layers will interact with each other the most if they were stacked in the AA 
arrangement (with all the carbon atoms on top of each other), leading to a more prominent band splitting, and thus 
the smallest gap possible from various stacking arrangements. For the AB stacking, however, in the a case, the 
difference with respect to the AA stacking is much less because the layers are only displaced in the x direction, while 
in the /3 case, layers are displaced both in the x as well as the y directions, resulting in more disalignment compared 
to the a case (c/. Fig. [1]). Therefore, in a alignment the layers will interact with each other more than those in the 
/? case, resulting in a band structure more distinct compared to a mono-layer, and, thus, a smaller gap, compared 
to the /3 case. We verified this hypothesis by actually performing the band structure calculations for 2-AGNR-ll in 
different stacking patterns and alignments, and the band gaps (in eV) obtained by using the TB/PPP-RHF method 
were 0.02/0.43 (AA), 0.04/0.92 (AB-a) and 0.13/1.12 (AB-/3). 

Next, we illustrate the variation in the band structure with the number of layers (n), by presenting the band 
structure near Ep of 3-AGNR-8-a (FigEli), 4-AGNR-8-a (FigEb), and lO-AGNR-8-a (FigEt). The band structure 
changes significantly with increasing n, leading not just to more bands, but also causing bands near Ep to become 
flatter. Furthermore, this flatness extends more and more into BZ as n increases, which shows up as a tremendous 
increase in the joint density of states not just due to a denser band structure, but also due to van Hove singularity like 
increase caused by parallel bands. This has obvious implications for the optical absorption spectra of such ribbons, 
which will be discussed further in Sec. IIVDI 

C. Effect of gate bias on the electronic structure 

The variation of band gaps of bilayer AGNRs, when electric {E^) is applied along z-direction is quite interesting. 
Ab initio DFT study on these systems by Sahu et al.r^ revealed that the band gaps increase with the increasing E^ 
for ribbons with band gaps below a critical value, and decrease for ribbons with band gaps above this critical value. 
Using nearest neighbor TB approach, coupled with perturbation theory, they showed that for the a-aligned ribbons. 
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Figure 7: Band structure of (a) 3-AGNR-8-a , (b) 4-AGNR-8-a, (c) lO-AGNR-8-a near Ef obtained using the PPP-RHF 
approach. 
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Figure 8: (Color online) Variation of energy gap with for (a) 2-AGNR-8-a and 2-AGNR-24-a (b) 3-AGNR-8-a and 3- 
AGNR-24-a, obtained using the PPP-RHF approach. 

the critical band gap Ec = (v^ — l)t±, so that for t± = 0.4, i €c ~ 0.25 eV. In Fig. [8^ we present the variation 
of band for 2-AGNR-8-a and 2-AGNR-24-Q! with E^, obtained using the PPP-RHF calculations. In the absence of 
Ez, the band gap of 2-AGNR-8-a is larger than the ec, hence the band gap decreases with E^ (Fig. [8^, black solid 
line), whereas, the band gap of 2-AGNR-24-a in the absence of E^ is less than hence the band gap increases with 
Ez (Fig. [5^, red lines). In Fig. [5)3 we present the variation of band gap for S-AGNR-S-a and 3-AGNR-24-Q; with 
Ez , and we observe a trend similar to the case of the bilayer GNRs for 3- AGNR-8-a, but for 3- AGNR-24-a the gap 
exhibits small oscillations around a uniform value. The trends on the variation of band gaps with E^, obtained from 
our PPP-RHF calculations are fully consistent with the ab initio DFT results F^^ii^ 

Finally, we examine the variations in the band structures of multilayer AGNRs caused by the gate bias. In the 
absence of an external electric field, the fundamental gap of all the ribbons is located at A; = 0. However, when the 
electric filed is applied, location of the fundamental gap shifts in the BZ, from fc = to finite values of k, similar to 
what was observed in the gated multi-layer graphene as welli^ The rightward shift of fundamental gap continues with 
the increasing E^ which is evident from the Fig. [HI , where we have presented the band structure of of 3-AGNR-24-a 
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Figure 9: PPP-RHF band structure of 3-AGNR-24-a near Ef, for several values of external electric field 
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IV. OPTICAL ABSORPTION SPECTRA 



Basic formalism 



The optical absorption spectrum for incident radiation polarized in x or y or z direction is computed in the form of 
the corresponding components of the imaginary part of the dielectric constant tensor, i.e., eii{uj) using the standard 
formula 



.(c^) = C 



TT / a 



Mk)\pMm' 



dk. 



(3) 



where a is the ID lattice constant, |c(fc)) {\v{k))) denotes conduction (valence) band state, pi denotes the momentum 
operator in the i-th Cartesian direction, w represents the angular frequency of the incident radiation, Ecv{k) — 
ec{k) — £v{k), with edk) (e„(fc)) being the conduction (valence) band eigenvalues of the Fock matrix. 7 is the line 
width, while C includes rest of the constants. As in our previous work^the momentum matrix elements {c{k)\pi\v{k)) 
were computed numerically using the formula proposed by Pedersen et al^ 



B. Sensitivity of the optical absorption to Hamiltonian parameters 



Before we present and discuss our results on the optical absorption in multilayer ANGRs, we examine their sensitivity 
to the hopping and the Coulomb parameters used in the PPP model Hamiltonian. We used two sets of Coulomb 
parameters: U — 6.0 eV and Kij = 2.0 {i ^ j) and Ki^i = 1 (our modified screened parameters'^) and U = 8.0 eV 
and Ki,j = 2.0 {i 7^ j) and Ki^i = 1 (original screened parameters of Chandross and Mazumdar— ). Similarly, for the 
hoppings we used two sets: commonly used values t ~ 2.7 eV ,t = 0.27 eV, t± — 0.4 eV, tj^ — 0.3 eV, and a set 
with larger values of hoppings, viz., t = 3.16,^' = 0.316, t± = 0.39 eV, and 4 = 0.315.'*'' In Figs. HVBl and ITVB] 
we present the results of these calculations for the ribbons 2-AGNR-II-q; and 2-AGNR-12-Q;, respectively. From the 
figures the following trends emerge: (a) For 2-AGNR-II-q;, the change of hoppings (red line vs black line) does not 
change the results significantly, but the increase in the value of ?7 (J7 = 6 vs. U = 8) leads to a prominent blue 
shift in the spectrum, (b) For 2-AGNR-12-a, however, the situation is different because change in hopping leads to 
more significant changes in the spectrum, as compared to the change in U . Therefore, the dependence of the optical 
absorption on the Hamiltonian parameters is itself dependent on the nature of the multilayer ribbon concerned. 
Because, at present no experimental information is available as to the optical absorption spectra of GNRs, thus, it 
is impossible to judge as to what are the correct set of parameters to be used in the Hamiltonian. Therefore, given 
our past experience with the modified screened parameters {U = 6.0 eV and = 2.0 {i ^ j) and Ki^i = 1), and the 
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Figure 10: (Color online) Optical absorption spectra of 2-AGNR-ll-a calculated using the PPP-RHF approach for the light 
polarized along the x axis with Coulomb parameters (a) U = 6.0 eV and Kij = 2.0 (i 7^ j) and Hi^i — 1 and (b) U = 8.0 eV and 
Kij = 2.0 (i 7^ J) and Hi^i = 1. In both the cases the black lines correspond to hopping parameters (in eV) t = 2.7, t = 0.27, 
t\_ — 0.4 eV, and — 0.3, the red line corresponds to t — 3.16,t — 0.316, t\_ — 0.39 eV, and — 0.315. 
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Figure 11: (Color online) Optical absorption spectra of 2-AGNR-12-Q computed using the PPP-RHF approach for the lif 
polarized along the x axis, with various sets of PPP parameters as explained in the caption of Fig. IIV Bl 
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in-plane hopping t — 2.7 eV for the mono-layer GNRs,— in the remainder of this work we utilize modified screened 
parameters and the hopping values t — 2.7 eV ,t = 0.27 eV, t± = 0.4 eV, = 0.3 eV. 



C. Effect of edge alignment on the absorption spectrum 

In order to investigate the sensitivity of optical absorption to the type of edge alignment, we present the optical 
absorption spectra of 2-AGNR-8 in a and /? alignments. In Fig. [12] a we present the optical absorption spectrum for 
light polarized along the periodicity direction (x axis, exxi^))i where Sm„ denotes a peak in the spectrum due to a 
transition from m-th valence band (counted from the top) to the n-th conduction band (counted from the bottom). 
For a-alignment (black line), the first peak of exx{^^) at 1.06 eV is Sn, while the one at 1.96 eV represents £22- For 
the /? alignment (red line) Sn peak is located at 1.30 eV, while the second peak (E22) at 1.67 eV. The separation 
between the first two peaks is larger in the case of a-alignment, as compared to that in the f3 alignment. This is a 
consequence of the observation discussed in Sec. IIIIBI that the energy separation between the adjacent bands near 
Fermi energy at fc = is larger in the a-alignment compared to /? -alignment. 

In Fig. 112b . the absorption spectra corresponding to the y-polarized photons (ej,y(a;)) is presented, and the first 
peak for the a alignment (black line) located at 2.61 eV represents S14, the second peak at 2.69 eV is S22, while the 
peak at 3.3 eV is S51. Significant differences in tyy{i^) are observed for the /3 alignment: it starts with two tiny peaks, 
with the first one (Sn) at 1.30 eV, and the second one at 1.67 eV (S22), which were absent in case of the a alignment. 
The intense peaks for p -alignment are located at 2.62 eV (S^*), which is coincidental with the same peak for the a 
case, and 3 eV (S41+E23) which is located in the similar energy range as the higher peaks of the a case. 

Finally, we examine the z polarized component of the absorption spectra [e^zi^)) presented in Fig. 112b . for 2- 
AGNR-8. In the low energy regime ezzii^) follows exx{^) very closely , but the magnitude of the peaks of ezz(a;) is 
reduced significantly compared to the peaks of exx{^)j because of the weak inter-layer coupling. For a-alignment, the 




Figure 12: (Color online) Optical absorption spectra of 2-AGNR-8 calculated using the PPP-RHF approach for the light 
polarized along: (a) x axis (b) y axis and (c) z axis. The black lines represent the a alignment case, while the red lines denote 
the P alignment. A line width of 0.05 eV was assumed throughout. 



peak at 1.06 eV is Y}^ and second peak at 1.96 eV is E^^, while for the /3 alignment, the corresponding peaks are at 
1.30 eV, and 1.67, respectively. One very important difference in the spectra of the two types of ahgnments is that 
the peak intensities in the lower energy region are much higher in the a-alignment, than in the /3-alignment. 

Therefore, the absorption spectra for each polarization direction bring out distinct features depending on the type 
of edge alignment, which can be used in the experiments on oriented samples to determine the type of alignment. 



D. Variation of optical absorption with the number of layers 

Recently, Wright et alJ^ using the TB calculations performed on 3p + 2 family of bilayer AGNRs, showed that 
the inter-layer coupling causes strong inter subband transitions, leading to tremendous enhancement in the optical 
conductivity in the low-energy region, thus, making them suitable candidates for opto-electronic applications. In 
Fig. [T3]we present the optical absorption spectra for 1-AGNR-ll, 2-AGNR-II-q;, and S-AGNR-ll-a, calculated using 
the PPP-RHF approach for the x-polarized hght, with the aim of understanding the influence of n on the optical 
absorption. In 1-AGNR-ll, the peak of e^x{uj) is located at 1.31 eV, while in 2-AGNR-ll, and 3-AGNR-ll 
it becomes more intense, and successively gets red shifted to 0.94 eV and 0.84 eV, respectively. A similar trend is 
observed also for other transition peaks as well, as well as for the peaks in Cyy^Lu) and ezzi^) (figures not shown). 
Although, the band gaps obtained using the PPP-RHF approach are much larger than that using the TB method, 
yet the tendency towards the enhancement of the optical response with the increasing value of n in the lower energy 
region is obvious from these plots. This behavior is fully consistent with the earlier discussion in Sec. IIII Al that the 
enhancement of the optical response in the lower energy region of multilayer AGNRs can be understood in terms of 
the redshift of the fundamental gap, and the increase in the joint density of states near Ep, with increasing n. 



E. Effect of gate bias on the optical absorption 

As discussed in Sec. IIIIBl a gate bias has profound effects on the electronic structure of multilayer AGNRs, and, 
therefore, it is of great interest to imderstand the effect of electric field on the optical properties of these systems. In 
the FigHHa and [HI b we present the absorption spectrum for the a;-polarized light (exx(w)) for 2-AGNR-12-a and the 
2-AGNR-24-a respectively, for different values of . Upon increasing the value of E^ , we observe the following trends 
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Figure 13: (Color online) Optical absorption spectra of 1-AGNR-ll-a, 2-AGNR-ll-a and 3-AGNR-ll-Q calculated using the 
PPP-RHF approach, for the light polarized along the x axis. A line width of 0.05 eV was assumed throughout. 




Figure 14: (Color online)PPP-RHF absorption spectra at several values of the static electric field E^, for the incident light 
polarized along the x axis for (a) 2-AGNR-12-Q, (b) 2-AGNR-24-Q. A line width of 0.05 eV assumed throughout. 



in the plot. In 2-AGNR-12-a, peak En, which corresponds to the fundamental gap, gets red shifted, without much 
change in the intensity at Ez = IV/nm. But, for Ez = 2 V/nm, besides exhibiting the redshift, the peak also becomes 
more intense. In case of 2-AGNR-24-a, whose gap is less than Cc, the En peak expectedly blue shifts with increasing 
Ez, however, it loses significant intensity as the field is increased from Ez = to Ez = 1 V/nm. However, upon further 
increasing the field to .E^ = 2 V/nm, no significant change in the intensity is observed. Thus, as far as the intensity 
of Ell peak with increasing gate bias is concerned, 2-AGNR-12-Q! and 2-AGNR-24-Q; behave quite differently. To 
elucidate this point further, in FiglTS]we present the variation of the heights of the the En peaks of 2-AGNR-12-Q; 
and 2-AGNR-24-Q;, as functions of the gate field, in the range < Ez < 2.5 V/nm. Recalling that the height of a 
peak is proportional to the oscillator strength of the corresponding transition, we observe the following trends: (a) 
for 2-AGNR-12-Q! the intensity varies slowly for < Ez < 1.0 V/nm, thereafter it increases monotonically, while (b) 
for 2-AGNR-24-a the intensity decreases rapidly for < Ez < 1.0 V/nm, and, thereafter it exhibits slow oscillations 
around a constant value. Thus, the two ribbons exhibit opposite behavior as far as the dependence of the oscillator 
strength of En peak on the bias field is concerned. This behavior is in sharp contrast to the experimenta l^ ^ i and 
theoretical^ results on bilayer graphene, which report a monotonic increase in the intensity of the first peak, with 
an increasing gate bias. This peculiar behavior of the peak intensity as a function of the gate bias can be tested in 
oriented samples of multilayer AGNRs, and can be used to distinguish between the widths of such ribbons. 



V. SUMMARY AND OUTLOOK 



In summary, we have used the PPP model based 7r-electron approach, incorporating long-range Coulomb inter- 
actions, to study the electronic structure and optical properties of multilayer AGNRs, at the mean-field Hartree 
Fock level. We considered two types of edge alignments namely a and /3, both in Bernal (AB) stacking pattern, 
and systematically studied the variation of band gaps and optical absorption spectra of multilayer AGNRs with the 
increasing number of layers. These calculations have demonstrated anisotropic optical response of multilayer AGNRs, 
with the calculated absorption spectra being crucially dependent on the polarization direction of the incident light. 
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Figure 15: Variation of peak strength of En transition of e^x with static electric field E^, computed using the PPP-RHF 
method, for: (a) 2-AGNR-12-Q (b) 2-AGNR-24-a. 



Furthermore, the optical response has been shown to depend upon the type of the edge aUgnment {a vs. /3) and the 
number of layers, leading to the possible experimental determination of the nature of the edge alignment, and the 
number of layers in the system, using optical probes. Our calculations also reveal that a gate bias along the inter-layer 
direction not only alters the peak positions in the absorption spectra due to the change in the band gaps, but also 
the peak intensities (oscillator strengths) in a nontrivial way. The variation of the peak intensities as a function of 
the gate bias is of the opposite nature for the narrow and the broad bilayer AGNRs. a behavior in sharp contrast to 
bilayer graphene. These predictions of ours can be tested in optical experiments performed on the oriented samples 
of multilayer AGNRs. 

As far as the future studies on multilayer AGNRs are concerned, it will be interesting to probe the influence of 
electron-correlation effects on various properties studied here. Furthermore, the nature of excitons in multilayer 
AGNRs should be investigated, so as to provide a complete description of the linear optical absorption spectra of 
multilayer AGNRs. Calculations along those directions are underway in our group, and results will be reported in 
future publications. 
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